use TableA5.dta,clear


***column 1***
mixed happy  direct_election2002 nomination2002   age age_sq male  ///
ethnic  mood i.province||Village_code:, cov(un) mle
outreg2 using tableA5.doc, nolabel bdec(3) tdec(3)replace coefastr se  

***column 2***                                                                                              
mixed happy direct_election2002 nomination2002   age age_sq male married  ///
ethnic ccp cadre_now cadre_past edu_year  health mood l_houincpc ///
net_wealth working_hour unemployed i.province ||Village_code:, cov(un) mle
outreg2 using tableA5.doc, nolabel bdec(3) tdec(3)append coefastr se  

        
//the following two columns restrict the obs to respondents who are heads of household or heads' spouses 

***column 3***  
mixed happy  direct_election2002 nomination2002   age age_sq male  ///
ethnic ccp cadre_now cadre_past edu_year child health mood l_houincpc ///
net_wealth working_hour unemployed  ccp_s cadre_now_s cadre_past_s sp_edu_new ///
health_s unemployed_s   i.province if H2_101<=2 ||Village_code:, cov(un) mle
outreg2 using tableA5.doc, nolabel bdec(3) tdec(3) append coefastr se 

***column 4***  
mixed happy  direct_election2002 nomination2002 dist_town2002 minority incpc98 vpop2002  disaster2002  age age_sq male  ///
ethnic ccp cadre_now cadre_past edu_year child health mood l_houincpc ///
net_wealth working_hour unemployed  ccp_s cadre_now_s cadre_past_s sp_edu_new ///
health_s unemployed_s muchabove_aveinc above_aveinc muchbelow_aveinc below_aveinc ///
better_living worse_living expbig_increase expsmall_increase exp_decrease ///
 i.province if H2_101<=2||Village_code: , cov(un) mle
outreg2 using tableA5.doc, nolabel bdec(3) tdec(3) append coefastr se  

